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Abstract 


Meteorological droughts are a recurring natural phenomenon that causes lack of 
precipitation. The severity of meteorological droughts is estimated by 
established algorithms known as drought indices. One such procedure, perhaps 
the simplest, is the Reconnaissance Drought Index, or RDI, which is based on 
the ratio between precipitation and potential evapotranspiration (PET) for a 
determined continuous period of months. In this study, the RDI is applied to 
three durations of meteorological drought at a weather station selected from 
each of the three geographic or climate zones in the state of San Luis Potosi, 
Mexico, which are: Villa de Arriaga (Potosino Plateau), Río Verde (Mean Zone), 
and Xilitla (Huasteca Region). The monthly rainfall records and average and 
minimum temperatures of each station cover more than 50 years. PET was 
estimated by four methods: (1) the Penman-Monteith formula, which is the 
reference method, (2) the Thornthwaite, (3) the Turc, and (4) the Hargreaves- 
Samani. The operating procedures for these criteria are detailed in appendices. 
The analysis of the results indicates that the RDIs estimated with the 
Hargreaves-Samani method are best for reproducing the results of the Penman- 
Monteith formula, in the three climatic locations processed. The Turc method 
also led to results similar to those of the reference. Therefore, it can be said that 
the RDI is a robust drought index, which practically does not depend on the 
method of estimating the PET. There is a noticeable difference in the operational 
procedures of the Penman-Monteith formula and the Hargreaves-Samani 
method. The latter is a practical solution that is worth mentioning. 


Keywords: Meteorological droughts, potential evapotranspiration, statistical 
tests, mean square error, mean bias error, types of meteorological drought 
(light, moderate, severe and extreme). 
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Resumen 


Las sequías meteorológicas son un fenómeno natural recurrente que origina una 
escasez de precipitación. La severidad de las sequías meteorológicas se estima 
a través de algoritmos establecidos, conocidos como índices de sequías. Uno de 
tales procedimientos, quizá el más simple, es el índice de reconocimiento de 
sequías O RDI (Reconnaissance Drought Index), que está basado en el cociente 
entre la precipitación y evapotranspiración potencial (ETP), ocurridas en un 
cierto lapso, seguido de meses. En este estudio se aplica el RDI en tres 
duraciones de sequía meteorológica, en cada una de las tres estaciones 
climatológicas seleccionadas de cada zona geográfica o climática del estado de 
San Luis Potosí, México, que fueron: Villa de Arriaga (Altiplano Potosino), Río 
Verde (Zona Media) y Xilitla (Región Huasteca). Los registros mensuales de 
precipitación y temperaturas media y mínima de cada estación abarcan más de 
50 años. La ETP se estimó con cuatro métodos: 1) la fórmula de Penman- 
Monteith, que es el criterio de referencia; los criterios de 2) Thornthwaite; 3) 
Turc, y 4) Hargreaves-Samani. Los procedimientos operativos de estos criterios 
se exponen en los apéndices (ver más adelante). El análisis de los resultados 
indica que los RDI estimados con el método de Hargreaves-Samani es el que 
mejor reproduce los resultados de la fórmula de Penman-Monteith en las tres 
localidades climáticas procesadas. También el método de Turc conduce a 
resultados similares a los de referencia y por ello se puede establecer que el RDI 
es un índice de sequías robusto, que prácticamente no depende del método de 
estimación de la ETP. Al haber una diferencia notable en los procedimientos 
operativos de la fórmula de Penman-Monteith y del método de Hargreaves- 
Samani, este último es una solución práctica muy importante. 


Palabras clave: sequías meteorológicas, evapotranspiración potencial, pruebas 
estadísticas, error cuadrático medio, error de sesgo medio, tipos de sequías 
meteorológicas (ligeras, moderadas, severas y extremas). 
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Introduction 


The meteorological drought is a regional natural phenomenon produced by 
climate variability, which causes a decrease in the normal precipitation in an 
area over a significant period of time. For that reason, it has adverse effects on 
nature and society. The severity of meteorological droughts is commonly 
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estimated based on drought indices, which vary in complexity, ranging from 
those using a single climatic variable such as the SPI (Standardized Precipitation 
Index) to those developing a water-soil balance such as the PDSI (Palmer 
Drought Severity Index). Hao and Singh (2015) have found that a single variable 
is not sufficient to characterize droughts, because they are recurring natural 
phenomena caused by multiple factors. Therefore, they have proposed 
multivariate drought indexes. 


The /atent variable approach allows the development of multivariate indices (Hao 
8: Singh, 2015), consisting of establishing new climatic variables by means of a 
difference, or quotient, of other variables with widespread physical significance, 
such as monthly precipitation (P) and potential evapotranspiration (PET). When 
P-PET difference and the SPI operational algorithm were used, the SPEI was 
developed (Vicente-Serrano, Beguería, 8 López-Moreno, 2010) and when the 
P/PET ratio was used, the RDI Reconnaíissance Drought Index was established, 
whose operational procedure is quite simple (Tsakiris 4 Vangelis, 2005; Tsakiris, 
Tigkas, Vangelis, 8 Pangalou, 2007; Vangelis, Tigkas, € Tsakiris, 2013). 


The objective of this study was to present, in detail, the process of calculating 
the annual RDI, with durations of 3, 6, and 12 months. This procedure is applied 
to monthly rainfall data and average and minimum temperatures in three 
localities in the state of San Luis Potosí, representative of its three geographical 
or climatic zones. For this comparison of the RDI, the PET was estimated with 
four methods, whose detailed description is presented in appendices. These are: 
(1) the Penman-Monteith formula, which was the reference method; (2) 
Thornthwaite; (3) Turc, and (4) Hargreaves-Samani. The three results are 
analyzed and conclusions are formulated. 


Methods and materials 


The RDIs: equations 


The reconnaissance drought index (RDI) is initially calculated as the quotient 
between the accumulated monthly precipitation and the respective potential 
evapotranspiration, in k months considered for each study year ¡ (Tsakiris Y 
Vangelis, 2005; Tsakiris et al., 2007; Vangelis et al., 2013; Campos-Aranda, 
2014): 
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In the previous equation, k is the duration of the meteorological drought studied, 
jthe month considered, and í a range from 1 to NA, which is the number of 
years of the processed records (> 30). Since the magnitudes of a, can be 
represented probabilistically by the log-normal distribution, the standardized 
RDI values are easily obtained with the equation: 


rDIi, =2 (2) 


Oy 
in which: 
y = In(a,,) (3) 


In Equation (2), yis the arithmetic mean and «J,the standard deviation of the 
valuesy;. The positive values of the RDIs indicate wet periods and the negative 


ones are meteorological droughts, with the following severity: light up to -1.00, 
moderate ranging from -1.00 to -1.50, severe ranging from -1.50 to -2.00, and 
lastly, extreme less than -2.00. The common durations of k are 3, 6, 9, and 12 
months, where the first three relate to the months with the highest percentage 
of precipitation and the fourth to the period from January to December. 
Durations less than one year may also correspond to the period of crop growth 
or times of high demand. Campos-Aranda (2014) present a comparison between 
the RDIst and the SPLI. 


Estimation of the reference PET 


Towards the end of the 1970s, the Food and Agriculture Organization of the 
United Nations (FAO) formulated guidelines for estimating water demands for 
crops (Doorenbos € Pruit, 1977). Advances in research and more accurate 
assessments of the use of water by crops show that the Penman method, 
suggested by the FAO, often overestimates the requirements, and that the 
alternative empirical criteria presented, in a variable way, closely represent 
reality (Allen, Pereira, Raes, € Smith, 1998). 
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In May 1990, FAO organized a panel of experts and researchers, in collaboration 
with the International Irrigation and Drainage Commission (ICID) and the World 
Meteorological Organization (WMO), to review the methods for estimating crop 
demands and establish their modifications. The panel recommended the 
adoption of the Penman-Monteith formula as the standard method for estimating 
reference evapotranspiration (£To) and advised on the estimation procedures 
for its various meteorological parameters (Allen et a/., 1998). 


Appendix 1 presents the Penman-Monteith formula, including the procedures for 
estimating its parameters based on meteorological data. Appendix 2 describes 
the FAO recommendations for an application with exclusively climatic data. The 
latter converts the Penman-Monteith formula into an applicable and valid 
worldwide method for calculating and comparing ETo. Allen et al. (1998) 
indicated that it is preferable to apply the Penman-Monteith formula even with 
the Appendix 2 approach, rather than use any other empirical method. 


Empirical methods for estimating PET 


In Appendix 3, Equations A.20 to A.31 represent the operational procedures for 
three empirical methods for estimating potential monthly evapotranspiration ( 
PET), which are: Thornthwaite, Turc, and Hargreaves-Samanli. 


Processed climate records 


The state of San Luis Potosí can be divided into three climatic regions, which 
are: Altiplano Potosino Plateau, Mean Zone, and the Huasteca Region. The first 
has a semi-arid climate, the second is temperate-dry, and the third is warm- 
humid. In each of these regions, the weather stations with more extensive 
records and with the least number of missing monthly rainfall data and average 
and minimum temperatures were searched. The following three were selected: 
Villa de Arriaga, Río Verde, and Xilitla. In each of them, the few missing rainfall 
data were considered equal to the mode, estimated based on fitting the mixed 
gamma distribution to all the available monthly values (Campos-Aranda, 
2005a). The missing average and minimum temperatures data were estimated 
with an interpolation procedure that took into account the trend observed in the 
month before and after the missing value. 
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The average monthly values of precipitation (mm) and average and minimum 
temperatures (% C) of each processed weather station are listed in Table 1, as 
well as their respective recording periods. At the Villa de Arriaga station during 
the period from 2010 to 2014, average monthly temperatures values were used, 
due to the fact that the records covered until 2009. Figure 1 shows the location 
of the three weather stations in the state of San Luis Potosí. 
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Figure 1. Geographical location of the three weather stations processed of the 
state of San Luis Potosí, Mexico. 


Also, Table 1 shows the monthly average values for solar radiation (Rs;), 
obtained from the maps by Hernández, Tejeda-Martínez and Reyes (1991), on 
pages 65 to 76, for the three locations of the selected and processed weather 
stations. Also shown are the magnitudes of the average monthly potential 
evapotranspiration (mm) (PET;), obtained through the equations presented in 
Appendices 1 to 3, applied on a monthly basis. 


Table 1. Monthly average values of the climatic elements indicated in the 
three selected and processed weather stations in the state of San Luis Potosí, 
Mexico. 


Description: Jan Feb | Mar Apr May Jun Jul Aug Sept | Oct Non Dec | Annua 
| 


Weather station: Villa de Arriaga (Longitude 01% 23" OG. Latitude 21% 54' N. Altitude 2170 masl. NA = 53 years) 
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Precipitation 13. 7.6 7.0 10.4 30.9 57.2 | 72.1 56.8 | 62.7 25.4 5.6 8.9 | 357.7 

(1962- 0 

2014) 

Precipitation 2.1 2.0 2.9 8.6 | 16.0 20.2 15.9 17.5 7.1 1.6 2.5 100.0 

percentage 3.6 

Average 13. 13. 15.9 19.2 20.9 20.8 19.8 19.5 18.8 16.6 14.6 13. 17.2 

temp. 0 9 4 

(1962- 

2009) 

Minimum 7.0 9.5 11.2 11.4 11.3 11.0 10.5 8.4 6.0 8.4 

temp. 4.3 4.9 4.8 

(1962- 

2009) 

Solar 400 | 380 505 650 650 550 520 420 380 450 430 390 480 

radiation 

(cal/cm?/day 

) 

PET: 84. 87. 126. 160. 177. 157. 151. 133. 114. 115. 97.8 88. 1495. 

Penman- 1 8 5 6 7 1 6 6 3 5 7 4 

Monteith 

PET: 37. 39. 102. 101. 44.7 39. 

Thornthwait 6 3 57.6 | 82.1 3 4 98.1 92.9 | 77.5 60.1 1 832.6 

e 

PET: Turc 82. 75. 113. 156. 162. 138. 129. 105. 104. 94.1 82. 1340. 
7 8 5 9 6 8 0 5 95.2 4 1 7 

PET: 87. 77. 121. 166. 180. 146. 139. 112. 110. 96.0 86. 1420. 

Hargreaves- 7 4 2 5 1 9 9 0 96.1 3 4 5 

Samani 

Weather Station: Río Verde (Longitude 99% 59' OG. Latitude 21% 56' N. Altitude 987 masl. NA = 54 years) 

Precipitation 12. 10. 9.4 32.7 36.5 88.7 | 88.3 71.7 103. 44.2 15.4 12. 526.2 

(1961- 2 8 4 9 

2014) 

Precipitation 1.8 6.2 6.9 | 16.9 16.8 13.6 19.7 8.4 2.9 100.0 

percentage 2.3 2.0 2.5 

Average 16. 18. 21.7 24.6 26.4 26.1 25.0 25.1 23.9 21.8 19.0 17. 22.1 

temp. 2 3 0 

(1961- 

2014) 

Minimum 10. 12.6 15.7 18.2 19.0 18.3 18.3 17.7 15.1 12.1 14.6 

temp. 8.6 10) 9.7 

(1961- 

2014) 

Solar 375 | 350 440 550 540 540 530 480 350 390 375 350 430 

radiation 

(cal/cm?/day 

) 

PET: 78. 87. 130. 157. 170. 160. 157. 148. 107. 110. 93.0 84. 1485. 

Penman- 4 1 6 7 0 5 6 1 5 3 7 4 

Monteith 

PET: 36. 47. 119. 154. 147. 135. 131. 107. 54.3 41. 1140. 

Thornthwait 6 5 84.3 3 2 0 4 3 1 82.6 0 6 

e 

PET: Turc 88. 8l. 115. 149. 150. 149. 145. 132. 104. 95.0 84. 1394. 
0 2 7 0 5 8 0 6 98.2 1 9 1 

PET: 90. 81. 124. 162. 171. 165. 163. 147. 101. 110. 95.5 86. 1500. 

Hargreaves- 8 6 4 1 9 0 1 8 4 5 9 9 

Samani 


Weather station: Xilitla (Longitude 98% 59” OG. Latitude 21% 23" N. Altitude 630 masl. NA = 50 years) 
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Precipitation 62. 65. 72.5 115. 175. 373. | 432. | 429. 566. 292. 101. 59. 2746. 

(1965- 6 3 3 5 9 2 9 al 5 5 0 2 

2014) 

Precipitation 2.6 4.2 6.4 | 13.6 15.8 15.7 | 20.6 10.7 3.7 100.0 

percentage 2.3 2.4 2.1 

Average 17. 18. 21.4 | 24.2 | 25.9 | 26.2 | 25.6 | 25.9 | 25.0 | 23.1 20.3 18. 22.7 

temp. 4 7 3 

(1965- 

2014) 

Minimum 12. 13. 15.8 18.5 20.6 | 21.2 | 20.9 | 20.8 | 20.2 18.1 15.6 13. 17.6 

temp. 6 4 4 

(1965- 

2014) 

Solar 350 | 350 | 400 480 500 460 490 450 310 375 360 300 400 

radiation 

(cal/cm?/día 

PET: 69. 78. 106. 127. 142. 132. 137. 132. 98.4 105. 87.4 | 74. 1292. 

Penman- 1 3 1 9 9 1 9 2 9 3 5 

Monteith 

PET: 42. 48. 113. 146. 148. 142. 142. 119. 62.6 | 47. 1187. 

Thornthwait 1 2 80.2 3 3 9 8 1 1 94.9 2 8 

e 

PET: Turc 85. 82. 105. 130. 139. 129. 136. 126. 103. 94.3 76. 1300. 
8 1 8 9 3 8 1 6 89.9 0 8 3 

PET: 87. 82. 112. 140. 157. 141. 152. 141. 110. 95.1 77. 1390. 

Hargreaves- 9 5 4 2 2 1 9 3 92.2 0 2 0 

Samani 


Quantification of differences with MSE and MBE 


In order to quantify the numerical differences between the annual values of the 
RDIst of duration k, due to the effect of the change in the method of estimating 
the PET: the following two statistical indicators were applied: (1) mean square 


error (MSE) and (2) medium bias error (MBE); whose expressions are (Vangelis 
et al., 2013): 


ECM; =|2 


ESM, = NA (Xrefi — Xest) (5) 


, , / 
yA (Xrefi —xestg) |" (4) 


In the above expressions, Xref; are the annual values of RDIst calculated with 
Equations (1) to (3), for a duration k in months, with the PET, estimated using 
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the Penman-Monteith formula (Appendices 1 and 2), which is the reference 
method, and xXest; are the same annual values as the RDIst but calculated base 


on PET, , estimated by each one of the three empirical methods presented 
(Appendix 3). Equations (4) and (5) were also applied to the annual PE7" values 
(k=12). 


Analysis of results 


Homogeneity of climatic records 


Annual values were obtained based on each completed record of precipitation 
and monthly average and minimum temperature. With these series, the 
statistical quality analysis of the record was performed, for which the following 
seven tests were applied, one general and six specific: (1) Von Neumann, which 
detects loss of randomness by unspecified deterministic components, (2) 
Anderson and (3) Sneyers, which determine persistence, (4) Kendall and (5) 
Spearman, which detect trends, (6) Bartlett test of variability, and (7) Cramer, 
to identify changes in the mean. In all tests, a level of significance (a) of 5% 
was used. The statistical tests cited can be found in WMO (1971) and Machiwal 
and Jha (2012). The results of these tests are shown in Table 2, where NH and 
H represent non-homogeneous and homogeneous series or records, 
respectively. 


Regarding the annual precipitation (P), the records from Villa de Arriaga and Río 
Verde showed persistence, which was also detected with the von Neumann test. 
When taking into account that persistence is a statistical component of the time 
series, analyses aimed at quantifying meteorological droughts can continue, 
since the three records show no trend or changes in the mean, that is, loss of 
homogeneity. 


Table 2. Results of statistical tests applied to annual records of precipitation 
(P) and average (Tt) and minimum (t) temperatures from the weather stations 


studied. 
ds Villa de Arriaga Río Verde Xilitla 
Statistical tests: 
P Tt t P Tt t Pp Tt t 
1. Von Neumann NH NH NH NH NH NH H NH NH 
2. Persistence (Anderson) NH NH NH NH NH NH H NH NH 
3. Persistence (Sneyers) NH NH NH NH NH NH H NH NH 
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4. Trend (Kendall) H NH? H H NH? | NH? H NH? H 


5. Trend (Spearman) H NH? H H NH? | NH? H NH? H 
6. Variability (Bartlett) H NH NH H H H H H H 
7. Change in the mean (Cramer) H NH H H NH NH H NH NH 


First-order linear correlation coef. 1 (ri) | 0.523 | 0.791 | 0.494 | 0.285 | 0.446 | 0.418 | 0.049 | 0.522 | 0.470 


The opposite occurred with annual average temperature (Tt) records, which 
were totally non—-homogeneous, since they presented persistence (associated 
with r,), an upward trend, and a change in the mean. In relation to the minimum 
annual temperature (t) records, the Rio Verde station was the least 
homogeneous, since it had an upward trend, and Villa de Arriaga was the most 
homogeneous, since it only had persistence. Xilitla was also non-homogeneous, 
having persistence and a change in the mean. 


Numerical results of RDI:.t 


Based on the percentages of monthly precipitation shown in Table 1, the three- 
month period with more rain from July to September and the six-month period 
from May to October were defined. At Villa de Arriaga, Río Verde, and Xilitla, for 
k = 3 months the sums of percentages varied by 52% and for k = 6 months 
they varied by 83%. Once these periods were defined, Equations 1 to 3 were 
applied to the precipitation data and to the estimates of the PET. Due to space 
limitations, only some of these results are shown in Tables 3 to 5. 


Table 3. Data, RDIst values, and types of annual meteorological droughts (75) 
calculated with the estimated potential evapotranspiration, according to the 
empirical criteria indicated, in the Villa de Arriaga weather station, San Luis 


Potosí 
p According to Penman-Monteith formula According to Hargreaves-Samani method 
Year PET k=3 k=6 k=12 PET k=3 k=6 k= 12 
(mm) months months months months months months 
(mm) (mm) 
RDIst TS | RDIs | TS | RDIst | TS RDIst | TS | RDls | TS | RDIs | TS 
196 1495. - SL - SL - SL | 1433. - SL - SL - SL 
2 268.4 8 0.110 0.16 0.21 6 0.11 0.17 0.22 
9 7 2 9 9 
196 1442. = SL = SL = SL | 1387. = SL = SL e SL 
3 279.1 0 0.444 0.75 0.08 2 0.47 0.78 0.10 
5 9 4 6 8 
196 1392. - - - 1370. - - - 
4 489.0 8 0.600 0.96 0.91 2 0.55 0.88 0.85 
0 7 3 9 5 


196 1376. SS - e 1364. - - Sa 
5 426.1 5 0.600 0.38 0.70 0 0.55 0.30 0.63 
9 4 1 6 2 
196 1293. e = = 1287. = = - 
6 631.7 4 0.571 1.26 1.47 1 0.49 1.16 1.39 
4 4 5 1 5 1 
196 1299. - e - 1275. e e - 
7 538.0 7 0.722 0.99 1.19 1 0.66 0.91 1.13 
1 5 8 3 7 
196 1368. - - 5 1347. e e - 
8 379.9 5 0.664 0.58 0.52 8 0.59 0.50 0.45 
4 0 5 0 9 
196 1450. - SL - SL E SL | 1405. e SL e SL - SL 
9 174.0 0 0.285 0.78 0.89 9 0.33 0.84 0.92 
6 if: 1 1 4 
197 1418. e SL - e SL | 1381. = SL = SL = SL 
0 269.7 5 0.078 0.04 0.11 9 0.13 0.00 0.15 
5 9 3 9 9 
197 1394. - - - 1372. - - - 
1 593.0 1 1.259 1.57 1.24 1 1.19 1.48 1.17 
7 1 1 9 7 
197 1417 - - - 1388. - e a 
2 347.0 3 0.248 0.60 0.30 2] 0.19 0.53 0.25 
3 8 4 3 7 
197 1418. = - = 1390. - a = 
3 328.5 9 0.647 0.53 0.21 9 0.58 0.45 0.16 
8 3 6 8 2 
197 1423. - - SL e Ss 1391. - SL e Ss - Ss 
4 156.9 4 0.021 0.97 1.04 M 7 0.03 1.04 M 1.08 M 
8 1 1 4 1 
197 1432. - = e SL | 1397. == = - SL 
5 280.5 3 0.329 0.18 0.06 9 0.26 0.11 0.11 
5 9 7 8 2 
197 1397. - a - 1362. a a - 
6 359.5 9 0.695 0.74 0.39 3 0.63 0.66 0.34 
1 1 8 9 8 
197 1395. - SL e ss = SS | 1353. - SL = ss a ss 
7 96.0 6 0.319 1.53 1.83 3 0.33 1.56 1.86 
5 7 8 J 0 
197 1363. - - SS 1325. E - = 
8 507.5 3 1.321 1.33 1.01 2 1.28 1.27 0.97 
3 6 8 3 4 
197 1419 - - SL - SL | 1372. S > SL = SL 
9 192.0 7 0.077 0.77 0.69 de 0.05 0.80 0.71 
0 5 1 5 9 
198 1393. e e SL - 1355. e - SL - 
0 375.5 4 0.064 0.06 0.47 5 0.03 0.10 0.43 
2 0 9 1 0 
198 1385. - SL e SL SE SL | 1340. - SL = SL = SL 
1 270.9 6 0.567 0.34 0.07 1 0.57 0.36 0.10 
3 2 7 7 0 
198 1438. - SL = SL SL | 1376. = SL = SL a SL 
2 292.5 5 0.443 0.01 0.00 0 0.44 0.02 0.01 
5 6 9 4 5 
198 1385. - SL e SL - Ss 1345. - SL - SL - Ss 
3 154.5 6 0.219 0.79 1.02 M 3 0.23 0.80 1.05 M 
1 1 2 4 0 
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198 1419 3 SL - Ss - Ss 1372. - SL - Ss - Ss 
4 148.0 9 0.017 1.05 M 1.13 M 2 0.05 1.09 M 1.15 M 
5 5 9 0 6 
198 1428. a SL - SL = SL | 1369. = SL = SL = SL 
5 167.0 0 0.765 0.79 0.94 2 0.76 0.78 0.94 
2 1 0 9 9 
198 1404. - SL a - 1371. a SL e - 
6 504.0 8 0.621 1.01 0.95 9 0.63 0.97 0.90 
2 3 8 8 4 
198 1377. - = 5 1321. e e - 
7 504.9 7 0.939 1.04 0.98 5 0.94 1.05 0.97 
9 9 7 6 0 
198 1407. - - a SL | 1343. - - e SL 
8 256.0 2 0.337 0.05 0.19 9 0.33 0.51 0.20 
4 4 8 0 
198 1469. e SL - SL - 1354. 3 E SL - 
9 359.5 2 0.052 0.15 0.30 1 0.01 0.07 0.35 
5 7 8 4 8 
199 1507. - - - 1384. - E - 
0 728.3 6 0.916 1.44 1.45 9 0.98 1.51 1.50 
0 6 9 8 7 
199 1028. 1564. - - - 1416 z - - 
1 5 4 1.572 2.04 1.97 6 1.64 2.13 2.04 
3 7 9 8 9 
199 1624. = - = 1458. - a a 
2 838.0 5 0.015 0.93 1.56 3 0.11 1.05 1.65 
4 7 5 1 6 
199 1716. - - - 1639. - - e 
3 523.0 8 0.696 0.75 0.67 7 0.63 0.69 0.66 
4 7 4 0 7 
199 1561. - = - 1448. = = - 
4 564.9 4 0.632 1.12 0.96 4 0.70 1.19 1.00 
3 8 6 1 5 
199 1530. - SL - SL - 1409. e SL - SL - 
5 392.8 0 0.103 0.19 0.38 2 0.04 0.11 0.44 
0 8 2 5 0 
199 1686. - ES 5 1506. ES Sn = 
6 666.0 9 1.263 1.13 1.11 8 1.37 1.26 1.21 
6 5 0 0 5 
199 1730. a SL 3 SS 1536. =S - = 
7 454.0 8 0.052 0.27 0.42 6 0.04 0.39 0.53 
2 4 3 3 8 
199 1748. - e - 1564. a - - 
8 418.0 9 0.181 0.38 0.26 4 0.28 0.50 0.36 
3 7 5 4 9 
199 1698. - SL a SL - SL | 1516. - SL - SL - SL 
9 249.0 0 0.126 0.25 0.55 1 0.06 0.15 0.44 
6 8 5 3 9 
200 1599. - SE - Ss - Ss | 1469. 5 SE = Ss 5 SE 
0 123.0 0 5.157 1.46 M 1.64 4 5.14 1.41 M 1.58 
4 8 8 1 2 
200 1564. - SL e Ss - SS | 1458. = SL = Ss = ss 
1 111.8 2 0.470 1.48 M 1.77 a 0.45 1.43 M 1.73 
4 3 2 8 0 
200 1582. a SL - Ss e Ss 1516. - SL e Ss - Ss 
2 139.3 7 0.956 1.43 M 1.42 M 8 0.95 1.41 M 1.42 M 
7 1 6 9 6 
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200 1638. e - e 1590. e - - SL 
3 337.0 0.463 0.35 0.01 3 0.45 0.34 0.02 
6 4 5 3 1 
200 1645. - = = 1532. = = - 
4 613.4 0.452 0.97 1.01 7 0.48 1.01 1.04 
8 9 1 7 8 
200 1667. - - SL - SL | 1548. e - SL a SL 
5 287.0 0.147 0.32 0.28 7 0.19 0.26 0.24 
3 7 5 2 6 
200 1637. - = 5 1531 e e - 
6 340.0 0.120 0.12 0.02 1 0.15 0.15 0.05 
1 9 2 8 8 
200 1630. e SL e SL e SL | 1531. e SL a SL e SL 
7 270.0 0.564 0.47 0.35 0 0.53 0.43 0.32 
3 3 7 9 9 
200 1602. - - - 1499. eS - - 
8 546.4 1.214 1.17 0.86 5 1.24 1.20 0.89 
7 7 7 7 0 
200 1518. - SL e SL - SL | 1468. - SL - SL - SL 
9 201.0 0.501 0.62 0.73 8 0,49 0.61 0.75 
7 1 7 8 5 
201 1498. - - SL - SL | 1420. - - SL - SL 
0 264.0 0.029 0.76 0.24 1 0.03 0.75 0.24 
3 8 6 3 1 
201 1498. - Ss - SE - SE | 1420. e Ss - SE he SE 
1 75.8 1.324 M 2.07 2.35 1 1.32 M 2.06 2.33 
6 J 3 5 8 
201 1498. - SL - ss e Ss 1420. - SL e ss - Ss 
2 134.9 0.669 1.61 1.38 M 1 0.66 1.59 1.36 M 
0 2 5 9 9 
201 1498. - Ss - ss - Ss 1420. - Ss e ss - Ss 
3 138.9 1.436 M 1.84 1.33 M 1 1.43 M 1.83 1.32 M 
9 3 6 8 0 
201 1498. - ss - Ss - Ss 1420. a ss - Ss 3 Ss 
4 165.0 1.518 1.28 M 1.04 M 1 1.51 1.27 M 1.03 M 
7 3 7 6 0 


Table 4. Data, RDIst values and types of annual meteorological droughts (TS) 
calculated with the estimated potential evapotranspiration, according to the 
indicated empirical criteria, in the Villa de Arriaga weather station, San Luis 


Potosí 
Te ¿ According to Thornthwaite criterion According to Turc criterion 
a A A k=3 k=56 k = 12 PET k=3 k=56 k=12 
( ( months months months months months months 
) ) (mm) (mm) 
RDlst | TS | RDIs | TS | RDIst | TS RDlst | TS | RDIs | TS | RDIst | TS 
196 | 17. a SL a SL = SL | 1354.2 = SL 5 SL 5 SL 
2 4 9.0 | 825.3 | 0.12 0.20 0.21 0.11 0.17 0.22 
2 3 5 2 5 8 
196 | 16 776. - SL - SL - SL | 1319.6 - SL - SL = SL 
3 5 8.3 0.38 0.70 0.04 0.48 0.79 0.12 
3 2 5 5 9 0 
196 | 16 == = - 1309.5 = - A 
4 2 9.0 | 758.4 | 0.64 1.10 0.94 0.53 0.85 0.83 
4 7 6 8 9 0 
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196 | 15. Sl - a 1298.7 - = e 
5 9 8.9 | 756.8 | 0.62 0.41 0.71 0.53 0.28 0.61 
3 6 6 7 5 4 
196 | 14. = = - 1217.5 e E - 
6 0 7.5 | 700.3 | 0.63 1.33 1.51 0.47 1.13 1.37 
0 5 5 6 7 9 
196 | 13. - e - 1206.2 SN - a 
7 7 6.6 | 693.0 | 0.76 1.08 1.26 0.65 0.88 1.12 
6 9 0 1 6 6 
196 | 15. e e - 1282.4 - - S 
8 5 8.3 | 742.8 | 0.70 0.63 0.55 0.57 0.48 0.44 
9 9 3 9 0 3 
196 | 16. - SL | 0.85 | SL | 0.89 | SL | 1330.8 a SL - SL - SL 
9 8 8.8 | 798.5 | 0.30 2 4 0.33 0.84 0.92 
1 4 0 3 
197 | 16. = SL - = SL | 1312.1 = SL = SL E SL 
0 2 8.5 | 772.4 | 0.08 0.07 0.09 0.14 0.02 0.16 
6 l 4 3 6 7 
197 | 16. - - - 1309.9 e - - 
1 1 8.8 | 760.1 | 1.31 1.65 1.26 1.17 1.46 1.15 
0 8 9 2 0 1 
197 | 16. - - - 1321.0 e - - 
2 4 8.9 | 774.1 | 0.28 0.66 0.32 0.18 0.51 0.24 
4 4 9 1 1 2 
197 | 16. a e E 1323.0 - = = 
3 5 9.0 | 778.0 | 0.66 0.55 0.22 0.57 0.43 0.14 
9 5 8 1 8 8 
197 | 16. - - SL e = 1325.1 e SL - Ss 1.08 Ss 
4 6 9.0 | 780.5 | 0.00 0.98 1.03 M 0.03 1.05 M 8 M 
2 1 1 9 0 
197 | 16. e e = SL | 1331.4 - - = SL 
5 6 9.0 | 780.2 | 0.42 0.27 0.04 0.25 0.10 0.12 
3 1 5 1 1 6 
197 | 15. - - - 1300.8 - - - 
6 9 8.3 | 751.8 | 0.76 0.84 0.43 0.62 0.64 0.32 
5 8 9 0 2 7 
197 | 15. e SL - Ss 1.79 | SS | 1288.2 = SL = ss | 1.86 | SS 
7 6 7.9 | 7474 | 0.25 1.47 M 0 0.35 1.57 1 
7 5 1 8 
197 | 14. E - - 1259.3 a - e 
8 9 7.3 | 727.8 | 1.36 1.42 1.07 1.27 1.24 0.95 
5 8 8 0 4 7 
197 | 16. - > SL = SL | 1301.9 - = SL = SL 
9 0 8.2 | 768.0 | 0.92 0.74 0.66 0.04 0.81 0.72 
1 1 4 4 2 
198 | 15. - - SL e 1287.1 - - SL - 
0 7 8.2 | 754.5 | 0.05 0.05 0.50 0.29 0.11 0.41 
3 4 Z 7 7 
198 | 15. E SL = SL > SL | 1276.0 = SL - SL = SL 
1 4 7.6 | 737.9 | 0.50 0.25 0.00 0.58 0.38 0.11 
3 3 9 9 8 3 
198 | 16. = SL - E 1308.4 = SL 5 SL = SL 
2 1 7.9 | 766.2 | 0.37 0.06 0.05 0.46 0.04 0.02 
1 7 7 0 1 7 
198 | 15. - SL - SL - SL | 1279.0 a SL - SL S Ss 
3 4 7.7 | 739.9 | 0.08 0.67 0.96 0.24 0.81 1.05 M 
9 6 6 8 6 5 
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198 | 16. S = Isi - s | 1307.6 = si | - S - S 

4 1 18.2 |762.0 | 0.04 0.96 1.08 | M 0.07 1110 M| 116 | M 
3 8 9 3 3 4 

198 | 16. = Isi = [sil -= |sLÍ|i3016 = sil = [sil - [st 
5 o |7.8 | 764.1 | 0.72 0.75 0.88 0.76 0.79 0.95 
4 4 9 6 8 4 

198 | 16. -= [sl - - | 1303.5 -= [sl - - 
6 2 |8.8 | 767.6 | 0.62 1.05 0.97 0.64 0.95 0.88 
7 2 7 4 9 8 

198 | 14. z = - | 12447 - - 
7 8 | 6.8 | 732.4 | 0.94 1.08 1.05 0.93 1.03 0.96 
4 4 9 3 4 8 

198 | 15. - - - [sL| 1269.2 - a - | st 
8 2 | 6.9 | 745.8 | 0.36 0.08 0.12 0.32 0.03 0.19 
4 8 3 7 9 9 

198 | 15. = Isi] - J|sL - | 1281.7 - - | sl - 
9 6 | 5.9 | 755.8 | 0.01 0.05 0.43 0.01 0.08 0.35 
4 6 0 2 4 2 

199 | 16. - - - | 1299.8 - - - 
0 2 | 6.1 | 790.6 | 0.92 1.47 1.55 0.98 1.50 1.50 
6 2 1 5 9 7 

199 | 16. - - - | 1327.5 y - a 
1 9 | 5.7 | 821.0 | 1.62 2.08 2.07 1.64 2.13 2.04 
1 5 2 5 7 8 

199 | 17. -= [sl E - | 1362.0 E - z 
2 9 |6.2 | 881.9 | 0.09 0.85 1.60 0.14 1.07 1.66 
2 5 3 2 3 3 

199 | 22. | 12. | 1603. = Isi] - [| sL[o21 |sLuÍ|1454.2 - - 2 
3 2 8 7 0.37 0.58 0 0.78 0.84 0.76 
9 8 6 5 7 

199 | 17. z - - | 1345.7 - 3 - 
4 6 |7.6 | 872.3 | 0.56 0.97 0.95 2 0.72 1.21 1.02 
4 7 3 2 6 5 

199 | 16. = Isi] - |sL - | 1306.3 = sil - |sL 5 
5 5 | 6.3 | 821.4 | 0.00 0.26 0.43 0.04 0.07 0.46 
2 0 9 4 9 8 

199 | 19. E E - | 1387.1 - = E 
6 2 | 6.9 1|981.5 | 1.12 0.83 1.03 1.40 1.31 1.24 
6 7 2 1 5 9 

199 -= [sl - - | 1429.4 - E - 
7 | 19. | 7.0 | 996.7 | 0.24 0.41 0.35 0.87 0.44 0.55 
9 2 6 4 3 9 

199 | 20. 1055. - [st - - | 14511 - - a 
8 7 183 0 0.11 0.05 0.12 0.34 0.56 0.39 
2 6 0 1 1 6 

199 | 19. = Isi] - [| sL[o053 |sLÍ| 14208 = |sil - [sil - [sl 
9 3 | 6.7 | 923.3 | 0.06 0.23 2 0.06 0.13 0.43 
0 4 1 2 4 

200 | 18. sE S - [|ss| 1389.3 - [|se| - Ss - |ss 
0 3 |7.7 | 854.6 | 5.14 1.38 | M | 1.59 5.13 1.40 | M | 1.57 
4 5 7 3 2 4 

200 | 18. = Isi] - S - [ss| 1380.1 = |sit| - 3 = |ss 
1 1 18.2 | 845.1 | 0.40 1.44 | M | 1.74 0.45 1.43 | M | 1.72 
7 6 0 6 0 2 

200 | 19. | 10. = Isi] - S Ñ S | 1429.5 = [sil - S - S 

2 6 7 | 912.0 | 0.92 1.42 | M| 149 | M 0.95 1.40 M| 141 | M 
0 2 6 2 2 4 
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200 | 21. 13. | 1019. e - - SL | 1486.5 - = E SL 
3 5 1 0 0.48 0.34 0.18 0.46 0.35 0.00 
8 6 6 1 9 4 
200 | 20. 10. - = - 1446.7 = - - 
4 2 1 927.9 | 0.55 1.14 0.98 0.48 1.01 1.04 
5 8 8 1 1 2 
200 | 20. 10. - - SL > SL | 1458.5 5 SL - SL 
5 6 3 951.2 | 0.23 0.18 0.34 0.19 0.25 0.24 
4 2 2 8 8 0 
200 | 20. 10. e - = SL | 1446.0 = = = 
6 2 3 924.2 | 0.22 0.28 0.00 0.15 0.15 0.05 
8 6 6 1 5 7 
200 | 20. 10. - SL > SL e SL | 1446.0 e SL - SL - SL 
7 2 5 925.3 | 0.48 0.35 0.39 0.53 0.43 0.32 
5 0 9 4 6 8 
200 | 19. - e - 1409.1 - - - 
8 3 9.8 | 901.2 | 1.28 1.27 0.84 1.24 1.19 0.89 
9 8 1 4 9 3 
200 | 18. 10. - SL - SL - SL | 1386.1 - SL - SL - SL 
9 4 2 861.4 | 0.49 0.60 0,77 0,49 0.61 0.75 
8 6 8 6 6 0 
201 17. - - SL - SL | 1348.4 - - SL - SL 
0 2 8.4 | 808.6 | 0.04 0.74 0.20 0.03 0.75 0.24 
7 5 8 1 5 9 
201 17. - 5 - SE = SE | 1348.4 - Ss - SE a SE 
1 2 8.4 | 808.6 | 1.32 M | 2.08 2.32 1.32 M | 2.05 2.33 
6 3 4 3 7 2 
201 17. - SL 2 ss - Ss 1348.4 e SL - ss - Ss 
2 2 8.4 | 808.6 | 0.66 1.60 1.32 M 0.66 1.59 1.37 M 
1 8 4 8 5 0 
201 17. - Ss - SE = Ss 1348.4 + Ss 5 ss - Ss 
3 2 8.4 | 808.6 | 1.44 M 1.85 1.29 M 1.43 M 1.83 1.32 M 
0 2 7 5 2 1 
201 17. E ss - Ss e Ss 1348.4 - ss - Ss - Ss 
4 2 8.4 | 808.6 | 1.52 1.27 M 1.00 M 1.51 1.27 M 1.03 M 
3 9 5 7 4 3 


Table 5. Data, RDIst values and types of annual meteorological droughts (TS) 
calculated with the estimated potential evapotranspiration according to the 
indicated empirical criteria, in the Xilitla weather station, San Luis Potosí 


p According to Penman-Monteith formula According to Turc method 
Year PET k=3 k=6 k = 12 PET k=3 k=6 k = 12 
(mm) months months months months months months 
(mm) (mm) 
RDIst | TS | RDls | TS | RDIst | TS RDIst | TS RDIst TS | RDIst | TS 
196 3010. 1248. - - - 1293. - - - 
5 6 6 0.75 0.61 0.60 6 0.69 0.517 0.52 
5 6 6 2 2 
196 2627. 1228 - SE - SL - 1272. - SE - SL - 
6 7 8 2.20 0.10 0.12 7 2.31 0.214 0.00 
9 6 3 1 9 
196 2912 1280 - - - 1290 - - - 
7 1 6 0.64 0.44 0.37 1 0.60 0.501 0.39 
0 6 0 0 1 


196 2302. 1268. ES - SL - SL | 1289. - - SL - SL 
8 7 2 0.41 0.28 0.53 1 0.34 0.345 0.61 
0 8 6 8 2 
196 3532 1285. e = e 1301. = = - 
9 0 4 1.44 1.10 1.13 0 1.50 1.172 1.18 
8 2 3 4 2 
197 2961 1244. - e - 1273. e e - 
0 0 1 0.22 0.77 0.55 4 0.18 0.758 0.51 
7 9 4 9 8 
197 3330. 1301. - a E 1300. - - - 
1 5 2 0.74 0.94 0.84 9 0.73 0.951 0.93 
9 5 7 0 1 
197 3202. 1238. - e 5 1292. - SL e E 
2 0 5 0.07 0.77 0.88 6 0.01 0.639 0.78 
2 1 8 1 9 
197 3166. 1215 - e - 1290. - - - 
3 0 6 0.27 1.05 0.91 l: 0.12 0.911 0.74 
6 4 7 1 9 
197 3044. 1247. - e - 1290. - - - 
4 9 6 0.86 0.71 0.65 9 0.75 0.579 0.57 
5 0 5 4 9 
197 3474. 1294. - - - 1294. - - - 
5 5 7 1.30 1.18 1.03 0 1.51 1.377 1.13 
4 3 8 0 5 
197 3297. 1167. E = Sa 1263. a a = 
6 2 0 1.07 1.01 1.24 1 0.94 0.828 1.01 
3 4 5 8 4 
197 1787. 1255. - ss E ss 2 SS | 1285. - SE - ss - ss 
7 2 0 1.91 1.60 1.51 5 2.00 1.656 1.68 
3 1 6 8 6 
197 3276. 1234. - = - 1287. = = - 
8 6 0 0.78 0.90 0.99 7 0.76 0.886 0.90 
1 9 5 1 4 
197 2364. 1326. E SL - SL - SL | 1297. - SL - SL e SL 
9 1 2 0.11 0.69 0.61 6 0.02 0.622 0.52 
1 2 1 5 7 
198 1896. 1363. S SL - ss = ss | 1314. - SL - SE = ss 
0 6 4 0.61 1.79 1.61 2 0.44 14.70 1.52 
7 0 0 4 4 Y 
198 3351 1307. - > a 1304. - - = 
1 7 4 0.18 0.62 0.85 6 0.24 0.713 0.94 
3 he 3 9 6 
198 1951 1345. - ss - ss - Ss 1310. - ss - SE = Ss 
2 0 8 1.81 1.77 1.44 M 0 1.66 1.752 1.39 M 
2 0 4 7 2 
198 3728. 1320. - e - 1296. - - - 
3 5 8 1.79 1.39 1.24 7 1.90 1.561 1.42 
0 1 2 9 8 
198 3758. 1257. S e = 1287. e e = 
4 7 0 1.47 1.72 1.47 5 1.40 1.718 1.49 
3 0 4 0 3 
198 2720. 1296. - SL - SL = 1294. = SL - SL = 
5 6 9 0.44 0.17 0.04 7 0.48 0.222 0.08 
9 9 5 2 4 
198 2552 1284. - Ss - SL e SL | 1299. a Ss - SL - SL 
6 2 9 1.43 M 0.24 0.17 E 1.44 M 0.299 0.20 
3 4 5 9 5 
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198 2388. 1257. = 0.07 | SL - SL | 1272. - - SL - SL 
7 8 9 0.31 6 0.35 9 0.31 0.097 0.40 
9 6 0 1 
198 2815 1253. - SL = - 1288. E SL e - 
8 8 0 0.16 0.22 0.32 9 0.19 0.123 0.25 
7 7 3 8 1 
198 2713 1333. e SL - SL - SL | 1303. ee SL e SL - 
9 6 10) 0.38 0.85 0.07 6 0.42 0.859 0.04 
6 0 6 3 4 
199 2639 1285. - - SL - SL | 1302. - - SL = SL 
0 0 6 0.29 0.23 0.04 9 0.23 0.291 0.07 
0 1 2 2 4 
199 3597. 1256. - - A 1301. e e - 
1 0 3 0.85 1.19 1.29 9 0.77 1.194 1.25 
3 6 9 9 7 
199 3175. 1191. - e - 1279. - - - 
2 5 0 0.32 0.58 1.01 9 0.23 0.451 0.79 
0 5 2 6 6 
199 3481 1240. - e - 1296. - - - 
3 9 3 0.55 1.01 1.21 0 0.47 0.922 1.13 
7 0 9 7 7 
199 2532 1271. - e - SL | 1305. - - SL - SL 
4 3 3 0.33 0.02 0.16 9 0.35 0.019 0.26 
3 9 4 1 0 
199 2617 1292. e S - SL | 1313. e ES e SL 
5 0 1 0.36 0.00 0.09 3 0.35 0.009 0.14 
2 2 6 6 3 
199 1918. 1305. - SL e Ss E Ss 1304. a SL e Ss - Ss 
6 0 1 0.69 1.08 M 1.38 M 4 0.71 1.107 M 1.44 M 
5 6 9 3: 6 
199 1999. 1273. - SE e Ss e Ss 1302. - SE - Ss - Ss 
7 5 9 2.12 1.26 M 1.12 M 9 2.15 1.344 M 1.26 M 
6 9 3 6 3 
199 2819 1314. - a - 1326. e e - 
8 5 2 0.56 0.15 0.13 8 0.57 0.201 0.13 
7 5 6 2 2 
199 2576. 1312 - = - sL | 1313 = A - SL 
9 1 0 0.72 0.18 0.22 4 0.71 0.198 0.21 
9 8 1 0 1 
200 1974. 1307. a ss - Ss - Ss 1317. = ss si Ss e Ss 
0 6 4 1.86 1.21 M 1.27 M 6 1.89 1.309 M 1.36 M 
0 4 9 8 5 
200 2751 1293. e SL - SL - 1313 = SL = SL - 
1 0 4 0.15 0.05 0.10 9 0.16 0.059 0.06 
8 3 1 1 9 
200 1880. 1314. E SL - Ss - Ss 1309. - SL - SL - ss 
2 7 4 0.90 1.04 M 1.49 M 7 0.88 0.996 1.54 
2 0 7 8 8 
200 2686. 1294. S E a 1305. - e z SL 
3 6 2 0.55 0.20 0.00 0 0.44 0.187 0.00 
7 4 3 8 4 
200 2155. 1283. E SL - SL = SL | 1301. - SL = Ss > SL 
4 8 1 0.75 0.96 0.84 5 0.76 1.041 M 0.93 
4 7 9 1 6 
200 2517 1332. - SL e SL e SL | 1316. - SL - SL - SL 
5 9 2 0.60 0.15 0.37 0 0.50 0.076 0.31 
1 2 5 7 8 
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200 1654. 1379 - Ss E SE - SE | 1322. Ss - SE - SE 
6 4 0 1.44 M 2.44 2.20 8 1.37 M 2.422 2.14 
3 8 6 8 0 
200 2948 1311 e = e 1305 = = - 
7 9 1 0.42 0.41 0.32 6 0.48 0.499 0.39 
0 0 6 6 4 
200 3475 1332 - e - 1310 a e - 
8 2 7 0.94 1.14 0.92 1 0.94 1.247 1.08 
3 3 2 0 3 
200 2239. 1393. e SL - SL - Ss 1325. e SL 5 SL = SL 
9 0 9 0.89 0.88 1.03 M 3 0.71 0.688 0.85 
2 7 0 4 1 
201 2880 1340 - > SL 5 1302 e E - 
0 7 4 0.71 0.01 0.14 2 0.80 0.102 0.30 
0 6 3 5 4 
201 1554. 1387 Ss SE A SE | 1313. Ss = SE = SE 
1 5 8 1.16 M 2.12 2.48 4 1.10 M 2.045 2.37 
9 2 3 7 6 
201 2202. 1376. - SL - SL - Ss 1315. - SL - SL - SL 
2 8 5 0.19 0.98 1.04 M 2 0.10 0.905 0.88 
3 7 6 8 8 
201 3764 1327 - e - 1304 - - - 
3 8 5 1.33 1.18 1.26 0 1.41 1.297 1.44 
3 6 0 0 6 
201 3100. 1356 - SL S = 1311. - SL e - 
4 4 1 0.45 0.46 0.39 7 0.41 0.530 0.58 
2 4 2 8 8 


Tables 3 and 4 show all the results of the Villa de Arriaga station, which is the 
one with the greatest variability, since its average annual precipitation and 
temperature ranged from 75.8 to 1 028.5 mm and from 13.7 to 21.5 *C. Table 
5 shows some of the results corresponding to the Xilitla station, which is the one 
with less dispersion, with average annual precipitation and temperature ranging 
from 1 554.5 to 3 764.8 mm and from 21.2 to 23.9 *C, respectively. 


Tables 3 to 5 use the following symbols for the severity or types of 
meteorological drought: light droughts (SL), moderate droughts (SM), severe 
droughts (SS), and extreme droughts (SE). The numerical results of the RDIst 
shown in Tables 3 to 5 allow a precise or detailed inspection and comparison of 
their annual values, observing a remarkable similarity both in their annual values 
and in the types of meteorological droughts they define, independently of the 
PET estimation method. The above will be numerically modified in Tables 6 and 
7, and can be seen in Figure 2 with the results in Table 3, relating to the Villa de 
Arriaga station, with k = 12 months. 


Table 6. Comparison of the MSE and the MBE between the annual PET and 
RDIst with the Penman-Monteith formula and their respective values as estimated 
with the three empirical methods cited, for the three indicated weather stations in 
the state of San Luis Potosí, Mexico. Note: minimum values of each comparison are 
shown in parenthesis. 


Station: Villa de Arriaga Station: Río Verde Station: Xilitla 
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Concept: 


MSE of PET 
annual de 
Thornthwaite 


MSE of PET 
annual de 
Turc 


166.8 


94.2 


(39.1) 


MSE of PET 
annual de 
Hargreaves- 
Samani 


(89.1) 


(24.5) 


103.1 


MBE of PET 
annual de 
Thornthwaite 


662.8 


344.8 


104.8 


MBE of PET 
annual de 
Turc 


154.7 


(-7.8) 


MBE of PET 
annual de 
Hargreaves- 
Samani 


(74.9) 


-97.5 


MSE of RDIst 
annuales de 
Thornthwaite 


0.165 


0.210 


0.136 


0.072 


0.100 


0.097 


0.086 


0.117 


0.112 


MSE of RDIst 
annuales de 
Turc 


0.062 


0.080 


0.062 


0.035 


0.044 


0.046 


0.082 


0.089 


0.108 


MSE of RDIst 
annuales de 
Hargreaves- 
Samani 


(0.050) 


(0.062) 


(0.050) 


(0.029) 


(0.036) 


(0.039) 


(0.070) 


(0.076) 


(0.098) 


MBE of RDIst 
annuales de 
Thornthwaite 


0.263-107 
e 


3.104-107 
7 


0.157-107 
” 


0.464-107 
7 


0.072:107 
” 


2.390-107 
7 


1.571-107 
) 


10.82-107 
7 


3.034-107 
7 


MBE of RDIst 
annuales de 
Turc 


1.468-107 
7 


2.525:107 
7 


1.338-107 
7 


3.024-107 
7 


1.352-107 
7 


4.136-107 
7 


4.908-107 
7 


3.648:-107 
7 


2.730:107 
) 


MBE of RDIst 
annuales de 
Hargreaves- 
Samani 


0.461-107 
7 


1.091:107 
a 


1.343-107 
7 


0.248-107 
) 


5.061-107 
7 


( 
1.209-107 


8.440:107 
7 


( 
0.071-107 


4.518-107 
7 


Table 7. Severity of the meteorological droughts obtained with the RDIst for 
the three durations (k) studied, in months, applying each one of the estimation 
PET criteria, in the three weather stations indicated, in the state of San Luis 
Potosí, Mexico. 


Types Penman-Monteith Thornthwaite Turc Hargreaves-Samani 
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hts 
(SMET 
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Weather Station: Villa de Arriaga (NA/2 = 26.5) 
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(RDIst)r-m 


(RDis)u 


(k = 12 meses) 


Figura 2. Comparison of the 53 RDIst values calculated with the PET based on 
Penman-Monteith (ordinate) and Hargreaves-Samani (abscissa), in the Villa 
de Arriaga weather station, San Luis Potosí 


Results of MSE and MBE 


Table 6 shows the numerical values of the MSE and the MBE. The comparison at 
the annual level of the PET estimates indicates that, in semi-arid and temperate- 
dry climates, the results of the Thornthwaite method is least similar to the 
Penman-Monteith formula, and the Hargreaves-Samani is more accurate. In 
the warm-humid climate, the results from the above two methods are nearly 
the same, and the Turc is most accurate. The minus sign in the MBE 
corresponding to the Turc and Hargreaves-Samani methods (last column of 
Table 6) indicates that these criteria overestimated the PET, with respect to that 
of the reference. The previous findings define the presentation of the results of 
the RDIs in Table 3, Table 4, and Table 5. 


The MSE corresponding to the annual RDls values for the three durations 
analyzed was greater with the Thornthwaite method and of similar order of 
magnitude with the other two criteria, but the Hargreaves-Samani method 
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always led to a lower value for the three climates studied. Regarding the MBE 
values obtained, in general they were low, of the same order of magnitude and 
modifying the results of the MBE, varying their sign according to the PET 
estimation method and the duration, k. 


Severity of meteorological droughts (SMET) 


Table 7 shows the estimates related to the number obtained from each type of 
SMET, for each of the three durations (k) and each weather stations processed. 
In general, the duration of three months showed the greatest dispersions in the 
SMET number, which in theory must be equal to half the number of years of 
records (NA), a value which is indicated for each weather station. The 
percentages quoted in Table 7 were calculated with the number of SMET found; 
therefore they add up to 100%. 


The Thornthwaite method resulted in the percentages of each type of SMET, 
which are more dissimilar than those obtained with the reference PET. This 
happened in the three weather stations, but was more pronounced in Río Verde. 


It can be stated that the percentages of each SMET that define the Turc and 
Hargreaves-Samani methods were quite similar to those obtained with the 
Penman-Monteith formula. This confirms the results in Table 6. 


The numerical values in Tables 6 and 7 allow us to conclude that there is no 
significant influence on the annual RDIst values, nor on the percentages of each 
type of SMET that they define, when the Hargreaves-Samani method is applied 
to any of the three weather stations processed. The Thornthwaite method is 
applicable only in the warm-humid climate of the Xilitla weather station. 


Conclusions 


The results of the application of RDIs to the three weather stations processed in 
the state of San Luis Potosí, belonging to different climates, indicate that there 
is no significant influence on the annual RDIst values or on the percentages of 
each type of meteorological drought that they detect, when using the empirical 
methods of Hargreaves-Samani and Turc to estimate monthly potential 
evapotranspiration (PET), in comparison with the results of the Penman- 
Monteith formula, taken as a reference. 
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This allows the RDIstto be established as a robust meteorological drought index, 
which is practically independent of the PET estimation method. 


The numerical calculation of the PET, according to the Penman-Monteith formula 
and the Hargreaves-Samani method, is remarkably different in terms of 
complexity. Therefore, the result in Table 6 indicating that the MSE is the lowest 
in the three climates studied, given such empirical criterion, is extremely 
important for its practical significance. 


Appendix 1: Penman-Monteith Formula 


Theoretical and operational equations 


H. L. Penman, in 1948, was the first to obtain an equation that combines the 
energy required to sustain evaporation and an empirical description of the 
diffusion mechanism by which energy is removed from the evaporation surface 
as water vapor (Shuttleworth, 1993). Penman's formula led to a new evaporation 
estimation criterion called the Combination Method. Several researchers 
modified the Penman formula to take into account the effects of the evolution of 
aerodynamic conditions on the growth of the crop, the former through resistance 
factors. The resistance of the surface (r,) on the water vapor flow in the stomata 


of the leaves and on the soil surface is distinguished from the aerodynamic 
resistance (r, ) that occurs due to the friction of the air flow over the vegetable 


surface. Although the exchange processes in the vegetation layer are much more 
complex, the measurements and calculations of latent heat flow, AET, have 
shown a high correlation, at least for a uniform grass surface. With such 
modifications the theoretical Penman-Monteith formula was obtained (Allen et 
al., 1998): 


_ A(Rn-G)+Pa' Cp (es-e)/Ta 
AET = A+y(1+r5/Ta) (A.1) 


where AET is the speed of evapotranspiration in megajoule per m? per day 
(MJ/m2/d), A is the slope at one point of the saturation vapor curve versus the 
temperature in kilopascal per *C (kPa/*C), Rn is the net solar radiation in 
MJ/m2/d, G is the flow of heat from the ground in MJ/m2/d, p,is the average 


density of air at constant pressure in kg/m”, c,is specific heat of the air at 
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constant pressure in MJ/kg/*C, (es - e) is the vapor pressure deficit of air in kPa, 
r, is the aerodynamic resistance in s/m, y is the psychrometric constant in 


kPa/*C, and r, the surface resistance in s/m. 


When considering a hypothetical vegetation surface that is 12 cm high, with a 
fixed surface resistance of 70 s/m, and an albedo of 0.23 in active growth that 
completely shades the ground and does not lack water, the following operational 
Penman-Monteith formula (Allen et a/., 1998) is obtained: 


_0.408-A4-(Rn—G)+y:[900/(Tt+273)]-uz"(es—e) (A.2) 


Elo A+ y:(1+0.34:uz) 


where ETo is the reference evapotranspiration in millimeters per day (mm/d) 
and the two new terms are Tt, which is the average air temperature at 2 meters 
high in *C, and uz, which is the average wind speed at 2 m high in m/s. The 
value r,= 70 s/m corresponds to a moderately dry soil surface resulting from 


frequent irrigation, approximately weekly. 


To get Equation A.2 from A.1, the depth of water in mm/d can be expressed in 
terms of energy received per unit area. This energy refers to the heat needed to 
evaporate the specified water depth, and is known as /latent heat of evaporation 
(A), which is a function of the water temperature (Ta) and is calculated with the 
following equation in MJ/kg (Allen et a/., 1998): 


2 = 2.501 — 0.002361-Ta(A.3) 


Since the value of A does not change much with Ta, Ta = 20 *C is used, and 
then A is approximately 2.45 MJ/kg, that is, 2.45 MJ are required to evaporate 
one kilogram of water or one liter. Then, a 1 mm depth of water is equivalent to 
2.45 MJ/m?, since 1 mm per m? is a cubic decimeter, that is, one liter. The first 
numerical coefficient in Equation A.2 converts the radiation, expressed in 
MJ/m2/d, to evaporation, in mm/d, and is equivalent to the inverse value of A 
(1/4 = 0.408). 


Equation (A.2) can be applied at intervals of one day, ten days, one month or 
even the total duration of crop growth or one year. To obtain ETo in mm/h, the 
numerator in the rectangular parenthesis is changed to 37 and all the variables 
are per hour rather than per day. For verification of the results, in humid tropical 
regions with a moderate average temperature (Tt=20 *C), ETo varies from 3 to 
5 mm/d; and with a hot climate, (7t >30%C) it ranges from 5 to 7 mm/d, these 
intervals increase by one unit in the arid zones (Allen et a/., 1998). In Mexico, 
applications of Equation A.2 have already been done by González-Camacho, 
Cervantes-Osornio, Ojeda-Bustamante and López-Cruz (2008), and Chávez- 
Ramírez et al. (2013). 
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Estimation of parameters A and y 


All the expressions presented below are from Allen et a/. (1998) and are used to 
estimate the potential evapotranspiration (PET;) in month j of each year ¡, which 


is required for the application of Equation 1. The slope (A in kPa/*C) in the vapor 
pressure curve of saturation at a point relative to the average air temperature 
(TE) in *C, is calculated with the expression: 


i 
17.27 'Tt; 
J 


4 089:|0.6108:exp - 
Ttl4237.3 


5 daa) 


Al = 
(rei+z37.3)' 


J 


The psychrometric constant (y in kPa/*C) is determined with the following 
expression: 


¡_ CpP _ 1.62862:10"3-P A 
Y EA; ax; QS) 


where c,= 1.013-1073 MJ/kg/*C and e = 0.622 are the quotients of the molecular 


weight of water vapor to that of air, A is estimated with Equation A.3 for the 
value of T+t; in *C, and P is the atmospheric pressure at the site in kPa. This is 


estimated with the equation: 
293 


.¿N 526 
P=1013: (=>) (A.6) 


in which, z is the altitude in meters above sea level. 
Estimation of radiations Rn and G 


The net radiation (Rn) is equivalent to the difference between the net short- 
wave incident solar radiation (Rns) and the net long-wave solar radiation that is 
emitted or released (Rn)), that is: 
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Rn = Rns — Rnl (A.7) 


The Rns is the difference between the incident radiation solar (Rs) and the 
reflected one, so it is estimated with the expression: 


Rnsj = (1 — a)-Rsj(A.8) 


where a is the albedo or coefficient of reflection of the vegetation cover, which 
is dimensionless. A value of 0.23 is adopted for the hypothetical reference grass. 
Rs; must be expressed in MJ/m?/d, therefore, the average monthly values, in 
cal/cm?/d, from the maps proposed by Almanza and López (1975), or by 
Hernández et al. (1991), must be multiplied by 0.041868 to obtain MJ/m/d. 


The long-wave energy emission rate is proportional to the absolute temperature 
of the surface raised to the fourth power. This relationship is known as the 
Stefan-Boltzmann's Law. Since water vapor, clouds, carbon dioxide, and dust 
absorb and emit long-wave radiation, their balance or net flow that leaves the 
earth's surface is estimated by correcting the Stefan-Boltzmann law for relative 
humidity and cloudiness, according to the following equation: 


Rnli = 0(Tt)- (0.34 - 0.14 [e)) : [135 : (52) - 0.35| (A.9) 


where, O = 4.903-107? MJ/K*/m/d is Stefan—Boltzmann's constant, Tt; is the 


average temperature of the month, in degrees Kelvin, equal to the degrees 


centigrade (*C) plus 273.16, e; is the current partial vapor pressure in kPa and 


Rso is the solar radiation on clear days or without cloudiness, in MJ/m/d. This 
is estimated with the expression: 


Rso = (0.75 + 2-10—5-z)-Rej(A.10) 


where, z is the altitude of the site and Re; is the what is known as extraterrestrial 
radiation, in MJ/m?/d. The quotient Rs/Rso must be less than one. The estimates 
of e;, Rej, and of the two missing terms in Equation A.2 (es and uz) are detailed 


in the following Appendix. 


By considering that the heat flow from the ground (G) is less than Rn, a very 
simple expression is used for its estimation, which considers the ground 
temperature to be similar to that of air; this ¡s: 
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a TEj+TEj-2 
SM 


Az (A.11) 


where, c.= 2.10 MJ/m3/“C is the caloric capacity of the ground, Ad is the interval 


in days, and Az is the soil depth affected, which for lapses of one month or more 
is considered equal to 2 meters. Based on these numerical values, Equation A.11 
for the first month, subsequent, and last month are: 


Gj=1 = 0.14 > (Tt;4, — Tt;) (A.12) 
Gj=12NA = 0.14 : (Tt; = Tt;_1) (A.14) 


In equation A.14, NA is the number of years in the climatic record processed. 


Appendix 2: Complementary climatic estimations 


Extraterrestrial radiation 


Rejy is the solar radiation at the top of the atmosphere in cal/cm?/d. It is tabulated 
monthly and is a function of the latitude of the site (q), in degrees. To avoid 
interpolation of Re;, 12 third-degree Newton polynomials were developed. Their 
formula is applicable at latitudes from 10 to 40 degrees north (Campos-Aranda, 
2005a): 


Re; = bo + bi(p — 10) + b2(p — 10)(p — 20) + ba(p — 10)(p — 20)(p — 30) (A.15) 


Los b, coefficients are as follows: 


Months | bo bi b> b3 Mes | bo bi b> b3 
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Jan 760 |-12|-0.075 1/600| Jul | 880 5 |-0.100 | -1/1200 
Feb |820| -9 |-0.100| 1/1200 | Aug |890 2 |-0.125| -1/1200 
Mar |875| -5 |-0.125| 1/1200 |Sept|880| -2.5 |-0.150| 1/1200 
Apr [|895| 0 |-0.125| -1/1200 | Oct |840| -8 |-0.075|-1/1200 
May |890| 4 |-0.100| -1/400 | Nov | 780 | -11.5|-0.025| -1/300 
Jun |875| 6 |-0.100 | -1/600 Dec 740 -12.5|-0.075| 1/1200 


The values of Re, estimated with Equation A.15 must be multiplied by 0.041868 
to obtain them in MJ/m?/d. 


Partial vapor pressure 


To estimate e; with Equation A.9, we remember that hotter air may contain more 


water vapor, whose maximum is the partial pressure of saturation vapor (es), in 
kilopascal (kPa), which is a function of temperature and is estimated with the 
following expression: 


17.27-Tt; 


es; = 0.6108- exp | | (A.16) 


(Ttj+237.3) 


When there is less amount of water vapor than the maximum, the partial 
pressure of water vapor is designated by e, and then the relative humidity (HR) 
in percentage is: 


HR = a 100 (A.17) 


If the air that contains any amount of water vapor equal to e is cooling, it reaches 
a point where e becomes es, and that temperature is called the dew point (+, ), 


which is estimated with Equation A.16. Solving for this gives us: 


alos 237.3|rm(ej/0.6108)] 
J 17.27-[Im(e¿/0.6108)| 


(A.18) 
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Having found that +, is very close to the average monthly minimum temperature 


(+,), a simple way of estimating the value of e; is established. This was verified 


almost three decades ago by Arteaga-Ramírez (1989), and more recently by 
Cervantes-Osornio, Arteaga-Ramírez, Vázquez-Peña, Ojeda-Bustamante and 
Quevedo-Nolasco (2013). 


Table A.1 shows the monthly average relative humidity data (Equation A.17), 
partial vapor pressure (e), and minimum temperature (£) for five meteorological 
observatories (SARH, 1982), which are surrounding the three weather stations 
that will be processed. Based on Equation A.18, using 6.108 instead of 0.6108, 
since e is in mbar, the corresponding dew point temperatures (+, ) were obtained. 
Then, the differences between t and +, were calculated for their inspection and 
to determine the corrections to the value t, since theoretically these differences 
should be close to zero. 


Table A.1. Monthly average values for several climatic elements in the five 
meteorological observatories indicated. *Average annual precipitation. 


Description: Ja | Fe | Mar | Apr | Ma | Jun | Jul | Aug | Sep | Oct | Nov | De | Annu 
n b y t Cc al 


Meteorological observatory: Saltillo (Coah.). PMA* = 269.4 mm. 


m 
N 


Relative humidity | 62 [| 59 | 54 | 54 | 58 | 62 | 65 | 68 | 72 | 70 | 64 62 


(Yo) 


Vapor pressure 91|9.|/9.00/13.|15. 17. 17. |17.5|16.|13.| 11. 9. | 13.3 

(mbar) 6 0 6 9 6 6 9 5 1 

Dew point (+,9C) | 5.6 | 6. | 5.4 | 10. | 13. |15.|15.|15.4| 14. | 11.| 90 | 5. | 11.2 
4 9 6 8 5 6 9 6 

Minimum temp. 5.2 | 6. | 8.7 | 12. | 14. | 16.|16.|16.2| 14. 11. 8.0 | 6. | 11.4 

(t *C) 8 7 7 4 5 5 6 3 

Differences £ —+, - |0.|33|18|1.1|06|1.0| 0.8 - - - O. 0.2 

(ec) 0.4 | 4 0.1|0.3| 1.0 | 7 

Meteorological observatory: San Luis Potosí (SLP). PMA = 315.4 mm. 

Relative humidity | 51 | 43 | 39 | 37 | 47 | 56 | 60 | 61 65 | 63 | 57 | 56 52 

(%) 

Vapor pressure 75|6.|7.3|7.9|11./|13.|13.|13.3|13. 12.|96| 8. | 10.4 

(mbar) 8 3 4 8 8 0 3 

Dew point (+.*C) |2.9| 1. | 2.5 | 3.6 | 8.8 | 11.|11.|11.2 11.|9.7/|6.4 | 4. 7.5 
5 3 8 8 3 

Minimum temp. 6.2| 7. | 9.9 | 11. | 13. | 14. |13.|13.5/| 13. |10.| 8.1 | 6. | 10.7 

(t *C) 4 9 4 2 5 2 7 5 

Differences t -—t, |3.3| 5.|7.4|8.3|4.6|2.9|1.7| 23 14|1.0|1.7 2. 3.2 

(ec) a 2 
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Meteorological observatory: Río Verde (SLP). PMA = 484.9 mm. 


Relative humidity | 73 | 70 | 64 | 64 | 66 | 69 | 73 | 72 | 77 | 76 | 76 [75| 71 
(Yo) 
Vapor pressure | 12.[13| 14. | 17. | 19. |21.[20.[21.4 | 20. |18.| 15. | 13 | 17.5 
(mbar) 0 0 3 09 E 2 E UK 946.7 
Dew point (+,9C) [10.[ 11 [| 12. | 15. | 17.| 18. 18. [18.6| 18. | 16. 13. | 11 | 15.4 
5 ea8s las Ba | EBe 4 po2 2 A 6 56 
Minimum temp. [9.2 | 10 | 12.| 16. | 17. | 19.| 18. [18.3 | 17.[15.| 12. | 9. | 14.7 
(E*0) 8|9|0|9|0 12 EA E 
Differencest=+. | - | - [(0.1[0.9[0.6|0.5|0.0|-0.3| - | - | - | - | -0.7 
(ec) 1.4 | 0. 0.9 |1.2 | 1.4 | 2. 
7 0 


Meteorological observatory: Aguascalientes (Ags.). PMA = 537.2 mm. 


Relative humidity | 57 [ 52 | 46 | 433 | 46 | 59 65 | 67 | 69 |64 | 59 | 61 57 
(Yo) 


Vapor pressure [8.8 | 8. | 9.4 | 10. | 12. | 14. | 15. [15.6 | 15. | 13.| 10. | 9. | 11.9 
(mbar) 9 0 A E E E E 

Dew point (+,9C) [5.1 | 5. [6.1 [7.1 | 9.9|12.|13.|13.6| 13. |10.| 7.7 | 6. | 9.5 
3 8 | 3 2 9 4 

Minimum temp. |4.6| 6. (8.2 | 11. | 13.| 14. | 14.[13.9 | 13. 10.| 7.2 | 5. | 10.2 
(t 90) 1 0 41 8|1 3 6 5 

Differences t=., | - [0.[|2.1[39|35|2.0|08| 0.3 |01| - | - | - | 07 
(ec) 0.5| 8 0.3| 0.5 | 0. 
9 


Meteorological observatory: Tampico (Tam.). PMA = 985.9 mm. 
Relative humidity | 81 | 81 | 80 | 82 | 81 | 82 | 80 | 80 81 | 79 | 79 | 80 80 
(%) 
Vapor pressure 17.|19 | 21. | 25. | 28.|30.| 30. | 30.4 | 29. | 25. | 21. | 18 | 24.9 


(mbar) 7 .3 1 5 5 5 3 7 9 5 9 
Dew point (+,*C) | 15. | 16 | 18. | 21. | 23. | 24. | 24. | 24.3 | 23. | 21. | 18. | 16 | 21.0 
6 9 4 4 2 4 3 9 7 7 6 
Minimum temp. 14. |15|17. | 20. | 22. | 23. | 23.|24.1|23.|21./| 18. | 15 | 20.1 
(t *C) 1 .7 6 8 9 9 8 1 2 3 1) 
Differences t —+, - - - - - - -= |-0.2 - - - - -0.9 
(90) 15|1.108|06/|03|0.5]|0.5 08 |0.5| 04 | 0, 
2 9 


The corrections in *C suggested in Saltillo are combined for the Villa de Arriaga 
station (September to February with zero and March to August with -1.00) and 
for Aguascalientes (July to February zero and March to June -2.00), proposing: 
September to February with zero and March to August with -1.50. Data from 
the San Luis Potosí observatory were not used since they are not considered to 
be fully reliable. For the observatory of Río Verde, the recommended corrections 
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are: March to August zero and September to February +1.20. Finally, in the 
Xilitla station, the correction suggested for the Tampico observatory, which is to 
add 0.60*C in all months, will be applied. 


Average wind speed 


Finally, regarding the estimation of the average wind speed (v) to be used in 
Equation A.2, the FAO has established four values according to the type of winds 
that occur in the region: (1) weak from 0.50 to 1.0 m/s, (2) moderate from 1 to 
3 m/s, (3) severe from 3 to 5 m/s and (4) strong > 5 m/s. It also points out that 
the average speed at 2 m high in more than 2 000 meteorological stations in the 
world is 2 m/s. 


Campos-Aranda (2005b) processed average monthly wind speed data (m/s) 
from 31 meteorological observatories that had this data, with the number of 
records ranging from 8 to 20 years. He found that the mode of 2.05 m/s verifies 
the observed global mean value and results in a sample and population median 
of the order of 2.3 m/s. This value is recommended for the monthly average. 


Solar radiation with temperature data 


When the Rs; maps of Almanza and López (1975), or Hernández et al. (1991) 
are not available, or one does not want to use monthly average values, an 
estimate based on the monthly difference between the maximum (7) and 
minimum temperature (t) of the air can be made, since such difference is related 
to the degree of cloudiness at the site. In general, on clear days or days without 
cloudiness, high temperatures are generated during the day and low at night, 
due to the fact that long-wave radiation is not absorbed or returned. The 
opposite occurs on cloudy days. The empirical operational equation is (Allen et 
al., 1998): 


Rs; = ka: Re; [T/ —tj (A.19) 


in which, Rs; is solar radiation, in MJ/m2/d, ka is an adjustment factor that ranges 


from 0.16 to 0.19, with units 1/%C, and Re; is extraterrestrial radiation in 

MJ/m2/d. The value ka = 0.16 is used for inland locations, where air masses are 

not influenced by the sea, and ka = 0.19 ¡s used in coastal areas where there is 
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such an affectation. Allen et al. (1998) also indicate how rs;data from a nearby 
meteorological station can be transported to the location under study. 


Appendix 3: Empirical estimation criteria for PET 


Thornthwaite Method 


This estimates the potential evapotranspiration (PET, ) of month j of year í based 
solely on the average temperature (T+;) in “C and the latitude of the site (LAT) 


in degrees. Its equation is (Mather, 1977; Campos-Aranda, 2005a, Xu, Singh, 
Chen, 8 Chen, 2008): 


Til 
10 Tt; 
IC; 


ETP| =16: ( ) "Fc; (A.20) 


in which, 1C¡ ¡is an annual heat index, equal to the sum of the 12 monthly indices, 
which are: 


The exponent ais a function of IC; with the following empirical equation: 
ar = 0.4924 + 1.792 : 107? + IC; — 7.71: 107? -IC2 + 6.75: 107 - IC? (A.22) 


Finally, Fc; is a monthly average corrective factor function of the latitude of the 
site and the number of days of the month (ndm). This formula is: 


Fc; = (2) . (==) (A.23) 
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where N is the maximum sunlight or maximum number of hours with average 
monthly sunshine. For its estimation in the Mexican Republic, Campos-Aranda 
(2005a) proposed the following empirical expression: 


N = A + B [sen (30 nm + 83.5)] (A.24) 


in which, nm is the number of the month, with 1 for January and 12 for 
December; A and B are constants, a function of the latitude of the site (LAT) in 
degrees, with the following expressions: 


> 
l 


12.09086 + 0.00266 - LAT (A.25) 


Uy 
ll 


0.2194 — 0.06988 - LAT (A.26) 
Turc Method 


In the early 1960s, L. Turc proposed the following equation to estimate the 
monthly and 10-day PET, which is a function of the average monthly 


temperature, T+t;, of the incident solar radiation Rs;, expressed in cal/cm?2/d, and 


of the average monthly relative humidity (Turc, 1961; Campos-Aranda, 2005a; 
Xu et al., 2008): 


ETP! =cC, 
J > Tt; Er 


E + 50): FC (A.27) 


The coefficient c; takes the following values: 0.40 for months of 30 or 31 days, 
0.37 for February, and 0.13 for a period of 10 days. The corrective factor is 
applied when the average monthly relative humidity (Hr;) is less than 50%, its 


expression is: 


FCi=1 + 11) (4.28) 
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Hargreaves-Samani Method 


Average daily potential evapotranspiration (ETP;),in millimeters, were proposed 
in the early 1980s, exclusively based on the average temperature (12) 


expressed in degrees Fahrenheit and the daily average incident solar radiation ( 
Rs;) expressed in millimeters of evaporated water depth. Its equation is 


(Hargreaves and Samani, 1982, Campos-Aranda, 2005a, Xu et al., 2008): 
ETP¡ = 0.0075 : Rs¡*Tt; (A.29) 


The incident solar radiation (Rs;) can be estimated with the Angstróm formula 


(Jáuregui-Ostos, 1978, Allen et al., 1998) when there is insolation or actual 
sunlight (n) data, or with the monthly maps available for the Mexican Republic 
(Almanza and López, 1975; Hernández et al., 1991), which are reported in 
cal/cm?/d. For the transformation of Rs, into evaporated water depth per day, 


the following formula is used: 


E 10:Rs; 
Rs] = Tp (4.30) 


where Ho, is the so-called latent heat of evaporation or energy needed in calories 


to evaporate 1 g or cm* of water. It is estimated with the following expression, 
with the average (T+;) monthly temperature in *C: 


Hvj = 595.9 — 0.55: Tt¡ (A.31) 


The other empirical formula by Hargreaves-Samani can be consulted in Campos- 
Aranda (2005a), which is a function of extraterrestrial solar radiation and 
monthly average and minimum temperatures. An application of this criterion is 
given in Campos-Aranda (2014). 
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